Origin of long-ranged attraction between like-charged particles at water-air interface 
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The nature of electrostatic interaction between like-charged particles at water-air interface is 
analyzed in this paper. We show that long-ranged electrostatic dipolar attraction between these 
objects generally exists. Our result provides a natural mechanism to explain the experimental 
observations of attraction between like-charge species trapped at water-air interface. We speculate 
that a similar (but presumably weaker) mechanism also exists for particles near water-solid interface. 
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When charged colloidal particles are dispersed into wa- 
ter, a fraction of the ionic functional groups on the sur- 
face dissociate. The net charge remaining on the particle 
surface is screened by the cloud of counter-ions in wa- 
ter, resulting in a screened Coulomb potential between 
particlespj 



U(r) 



-(r-2a)/X D 



e w (l + a/X D ) 2 



where r is the sphere's center-to-center separation, q is 
the effective charge carried by each particle with radius 
a(<< Xd), £w is the dielectric constant of water and 
Ad is the Debye-Hiickel screening length. This screened 
Coulomb repulsion together with short-ranged Van der 
Waals attraction usually determines the phase diagram 
of charged colloidal systems [3, [3]. 

In recent years, however, there has been increas- 
ing number of experiments reporting that under spe- 
cial circumstances, attractive interactions rather than the 
screened-Coulomb repulsion exist between like-charged 
species ranging from simple colloidal particles [3, 0, IE 
a, lZ| to complex cytoskeletal filamentous actin[8( and 
DNA0|. Because of its fundamental interest and im- 
portant implications in colloid science and biology, the 
problem has been under intensive theoretical study floL 
ll El EE El EE El for many years. 

In this paper, we shall study electrostatic interac- 
tion between charged particles trapped at water-air in- 
terface. It was discovered that colloidal particles that 
repel each other in deep water ma y a ttract each other 
when they are close to water-air 17, 18, IE 2ll or water- 
solid |, II 0, 0, interfaces. Moreover, whereas like- 
charge attraction between particles deep inside water 
usually occurs at distance scale r < A £> [E IE Ell j attrac- 
tion between colloidal particles near interface is found to 
appear at length scale r >> Xd 2l|. The long-range- ness 
of the interaction suggests that the origin of attraction is 
probably independent of details like nature of the screen- 
ing charge or chemical environment which largely deter- 
mines the interaction at length scale < Xrj. Indeed, it 
has been proposed [IE IS] that because one of the phases 
forming the interface (air or solid) is made of a non- 
polar substance that cannot sustain charge, the countcr- 




FIG. 1: A charged particle trapped in a water-air interface. 



ions in water are distributed asymmetrically around the 
particle, producing an effective dipole moment for each 
particle (Fig. [1]). For a spherical particle with uniform 
charge distribution on the particle surface, the resulting 
dipole moment points downward perpendicular to the in- 
terface, leading to long-ranged dipolar (~ 1/r 3 ) repulsion 
between particles. Following this line of argument, it was 
further suggested that charge inhomogeneity on the sur- 
face may lead to formation of an effective dipole moment 
not perpendicular to the interface, and ma y le ad to long- 
range dipolar attraction between particles [2 lj. 

Whereas this argument is physically suggestive, it has 
not been proven with rigor that the mechanism works. 
An electric dipole inside water is completely screened as 
monopole charges, and the interaction between dipolcs 
becomes short-ranged. It is not clear that the 1/r 3 
long-range interaction will survive for an effective elec- 
tric dipole trapped at the surface of water, or whether 
strong correction to the dipolar interaction exists. The 
purpose of our paper is to clarify this problem and show 
that the dipole picture is valid in the limit when A^> is 
smaller than any other length scales in the system and 
dipolar attraction between particles located at water-air 
interface appears under very general conditions that are 
independent of microscopic details of the system. 

We start by considering a particle of arbitrary shape 
trapped at a water-air interface (Fig. Q]). Charges from 
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the particle surface in contact with water dissociate leav- 
ing a surface layer of charge <j(x) at the particle-water 
interface. The dielectric constant of the particle is e p . 
Therefore the electric potential arising from nonzero cr{x) 
is obtained by solving the Poisson equation 



V 2 $(^ = 4ttV • P, 



at region I (air+particle) where 
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is the electric polarization inside the particle and 2?M = 
_ V$W(i = I, II) is the electric field at region i. The 
electric potential at region II (water) is obtained by solv- 
ing the linear Poisson-Boltzmann equation 
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subject to the boundary condition 

[e p E (I \x) - e w E {II) {x)] ■ n(x) = 4tt<j(x), (Id) 

where x runs over the surface S separating region I and 
II and n(x) is a unit vector normal to S (hereafter we will 
specify n pointing from regions II to I). Electric polariza- 
tion is introduced explicitly at region I for convenience as 
we shall see later. It is not introduced in water explicitly. 

Although exact solution to the above electrostatic 
problem is not available in general, the nature of the 
solution can be seen by using the Green's Theorem 
which states that the potential $^ / - ) at region I can be 
expressed in the form 
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where R ~ \x — x'\ and da' runs over the surface sepa- 
rating the two regions. The surface integral can be inter- 
preted as the potential due to a surface charge density 

a i = — 2^ e P^7 



and dipole density pi = j-Q^ri 24|. 



A similar theorem exists at region II with the 
Green's function 1/R for the Poisson equation in free 
space replaced by the corresponding Green's function 
^ e -R/\ D y i jj> f or t ne Poisson-Boltzmann equation, i.e. 
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Notice there is no "space charge" at region II. We note 
that the Green's theorem does not provide a full solu- 
tion to the boundary-value problem, but only an integral 



is an over-specification of the problem. Nevertheless, the 
theorem dictates the form of the solution. The resulting 
potential of any charge distribution a(x) on the surface 
separating the two regions can be viewed as coming from 
effective charge and dipole (normal to the surface) dis- 
tributions on the same surface plus induced electric po- 
larization inside the particles in free space. Interaction 
between particles are determined if we can find these ef- 
fective charge and dipole distributions. 

Eq. (fTc|) can be solved easily in the limit when \o is 
much less than the two other length scales in the prob- 
lem, the length scales where a{x) varies and the surface 
curvature of the interface. In this case, we can assume 
that the surface S is locally flat and cr(x) is locally con- 
stant. The solution to Eq. (fTc|) at region II not far from 
the interface is approximately 

$( // )(x', 2 /',z'))-*(^,y',0)e- |2 ' l/A -, 

where (x',y',0) denotes a point on the interface S and 
z' = zh' is a vector pointing away from the interface at 
point [x',y'). Correspondingly, 

^1 „ _ J_$(") 
dri Xrj 

on the interface. Putting this back into Eq. ^j, integrat- 
ing by part and making use of the boundary condition 
(lld[) . we obtain at region I 



(X) =-J v d 3 X 'P(x') ■ V'(i) + i § s {^ s {x')£r{\) 
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where = $( // ) = $g on the interface. <E>s can be 
determined easily if we impose the physical requirement 
that — > everywhere in region I in the metallic 

limit X D — > o[Il]. In this limit $s(f ) = Ait\ D (j{x')/e w 
and 



^\x) -> - / d 3 x'P(x'yv(l)+(( 



statement since arbitrary specification of both $ and §^ 4% f ^x E 1 (x) ■ E 2 (x) + \ D 2 &i(x)$ 2 (x) 



(4) 

corresponding to the potential from a volume elec- 
tric polarization P and surface dipole density p{x') = 
\D&(x')h(x')/e w . More generally, a systematic expan- 
sion of $5 in powers of A^ can be set up, and Eq. (|4]) 
gives the leading order result. Notice that Eq. ((3]) should 
be solved together with Eq. (|lb[) to determine P and 
<&( 7 ). Although analytical solution cannot be obtained 
in general, the dipolar nature of the electric potential of 
charge particles at the air- water interface is apparent. 

The interaction between two particles with 
surface charge distributions ui(x') and cr 2 (x') 
is given by the cross-term in the total electro- 
static energy, W12 = 37 J (7) d 3 xe p E 1 (x) • E 2 (x) + 

where 
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the first term is the contribution from region I and the 
second term from region II. The first term dominates in 
the limit — ► and 
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dvi / dv 2 (P 1 -V 1 )(P 2 -V 2 )G ( - I \x 1 ,x 2 ), (5) 
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where J v dvP.. includes both the volume and surface con- 
tributions to the electric dipole and 
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Xl\\X - x 2 \ 
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The first term in Eq. ^ gives rise to usually dipole- 
dipole interaction in free space whereas the second term is 
a correction term coming from the water-air interface. It 
is not difficult to show that n(x)-V(|x — xi||x — x^) -1 = 
for flat surface and the second term is generally smaller 
than the first term. The naive dipole picture [2ll I22I I23]] 
for interaction between colloids on water-air interface is 
recovered if we neglect the volume contribution P and 
the surface correction term in G(x*i,x 2 ). 

We next discuss corrections from finite value of A^. 
These corrections are generally small when the particles 
are far away from each other but may become large and 
dominant when they are within distance ~ Xjj . The most 
important contribution is perhaps the region II contribu- 
tion which is given by a formula similar to W± 2 , 
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where Oi is the surface charge density for particle i and pi 
is the corresponding surface dipole moment as discussed 
before. Notice that P does not contribute in this region. 
G 11 is the Green's function for the Poisson-Boltzmann 
equation, 



G( /7 )(x l7 x 2 ) = 
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The surface charges a[s with the first term in G^ 11 ^ give 
rise to the celebrated screened Coulomb interaction be- 
tween particles. The second term in G^ 11 ^ is a corre- 
sponding interface correction. 

Further corrections to effective dipole and surface 
charge distribution exist. These corrections reflect finite 
surface curvature or spatial dependence of a{x) and are 
important when the length scale is comparable to the 
screening length A^> . In this case the counter-ions in the 
fluid cannot adjust themselves to these rapid changes 
fully, resulting in formation of small patches of local 
charges or correction to magnitude of dipole moment at 
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effective surface charge density 
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uniform surface charge density g=g d 
curve surface z = h cos(27lx/a) re gion I 



region II 



effective dipole distribution 
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effective surface charge density 



(b) 

FIG. 2: Effective dipole distribution with dipole density 
pi — -^-&sn and effective surface charge with charge den- 



sity a i 



for two different situations, (a) A non- 
r) on a flat sur- 
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uniform charge distribution o~{x) — aod(ao 
face, all the charges appear inside a disk of radius ao, where 
r is the radial coordinate on the surface and ao = 5A_d. (b) 
An uniform charge distribution a(x) — <ro on a curve surface 
z = h cos(2nx / a), h — Xd and a — IOAd. 



length scale ~ \d- Notice that because of the long- 
rangcness of the Coulomb potential, non-local correc- 
tions to $s also arises and are important when surface 
curvature is nonzero. To illustrate we show in Fig. [2] 
numerical solutions of $W and dQ^^/dn for two dif- 
ferent situations: (i) a non-uniform charge distribution 
<t(x) = <7o0( a o — r ) on a na t surface, which is a two di- 
mensional disk (Fig. rjl^a)), and (ii) an uniform charge 
distribution a{x) = <7o on a curve surface S given by 
z(x,y) — hcos(2nx/a), which is stripe like (Fig. [2fb)). 
Water occupies the lower half of the figure (region II). 
We see that both surface curvature and non-uniform sur- 
face charge distribution give rise to corrections to $s 
and d&^/dn. For the first situation (Fig. Ufa)), the 



correction to effective dipole is Spi 



'-ao and the 



correction to effective surface charge 8a i ~ 777 t 2 - 00. For 



the second situation (Fig. [2Kb)), Spi ~ 2ir 



Ad ii 
2 z\o cr 



and 
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Summarizing, we study in this paper the nature of elec- 
trostatic interaction between charged particle trapped at 
the water-air interface. We find that besides screened 
Coulomb interaction, long-range dipolar interaction gen- 
erally exists between particles trapped at water-air inter- 
face, independent of microscopic details of the system. 
The interaction between particles can become attractive 
easily in the presence of non-uniform charge distribution 
on particle surface 21| , or for uniform charge distribution 
but when the shape of the particle is asymmetric (absence 
of inversion symmetry along one axis parallel to water 
surface) . In both cases, a nonzero average planar electric 
dipole moment may be induced, and leads to dipolar at- 
traction between particles. We note that a similar (but 
weaker) mechanism presumably exists for particles near 
a water-solid interface, if the distance between particle 
and the interface is less than the screening length \r> . In 
this case, the electric charge on the particle may not be 
totally screened and some of the electric field can leak 
out to the non-polar solid medium, leading to dipolar 
interaction between particles. We shall investigate this 
situation in more details in a separate paper. 

Lastly, we emphasize that the dipolar interaction we 
studied in this paper is limited to particle trapped at, 
or close to interfaces where the interaction is long-ranged 
in nature. It probably cannot provide a full explanation 
for like-charge attraction between particles deep in water, 
where interaction occurs mainly at length scale < Xd and 
many other factors contribute. We note also that there 
exists other mechanism which may induce long-ranged 
interaction between particles trapped at water-air inter- 
action, for example, van der Waals interacti onal and 
electric-field induced capillary attractions [3, lief . Our 
work provides a solid justification for the existence of 
the electric dipolar interaction. Which mechanism dom- 
inates depends on actual experimental conditions which 
probably differ for different experiments. 

We thank Dr. Wei Chen and Prof. Penger Tong for 
many helpful discussions and comments. This work is 
supported by HKUGC through grant CA05/06.SC04. 
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